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ABSTRACT 

We extract from the Sloan Digital Sky Survey a sample of 347 systems involving early type 
galaxies separated by less than 30 kpc, in projection, and 500 km/s in radial velocity. These 
close pairs are likely progenitors of dry mergers. The (optical) spectra is used to determine 
how the interaction affects the star formation history and nuclear activity of the galaxies. The 
emission lines (or lack thereof) are used to classify the sample into AGN, star forming or 
quiescent. Increased AGN activity and reduced star formation in early-type pairs that already 
appear to be interacting indicate that the merging process changes the nature of nebular ac- 
tivity, a finding that is also supported by an increase in AGN luminosity with decreasing pair 
separation. Recent star formation is studied on the absorption line spectra, both through prin- 
cipal component analysis as well as via a comparison of the spectra with composite stellar 
population models. We find that the level of recent star formation in close pairs is raised rela- 
tive to a control sample of early-type galaxies. This excess of residual star formation is found 
throughout the sample of close pairs and does not correlate with pair separation or with visual 
signs of interaction. Our findings are consistent with a scenario whereby the first stage of the 
encounter (involving the outer parts of the halos) trigger residual star formation, followed by 
a more efficient inflow towards the centre - switching to an AGN phase - after which the 
systems are quiescent. 

Key words: galaxies: elliptical and lenticular, cD - galaxies: evolution - galaxies: formation 
- galaxies: stellar content. 



1 INTRODUCTION 

It is well known that ear ly type galaxies are dominated by old stel- 
lar po pulations (see e.g. iKodama & Ar imoto 1997; Stanfor cfet al. I 
1998). It has also been shown through studies of NUV photome- 
try from GALEX dYi et al. II200I iKavirai et allbOOTh as well as 
statistical disections of the optical sp ectra dFerreras et al. I [20061 : 
iRogers et al. 1200 7; Nolan et al. 2007) that a large fraction of early 
types have undergone small amounts of recent star formation. How- 
ever there is less certainty as to the cause of this recent star for- 
mation. One possible scenario involve s minor mergers from small 
blobs of gas surrounding the galaxy (Kavirai 2008; Kavirai et al. I 
2009). In that case, the young stellar mass content will be roughly 
independent of galaxy mass, which implies recent star formation 
will be more readily detectable in lower mass early-types, as ob- 
served. This scenario will also result in enhanced recent star for- 
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mation within close pairs, where these small pockets of gas could 
be disrupted during the encounter. In this paper we study differ- 
ences in the spectroscopic data of close pairs involving only early- 
type galaxies with the aim of understanding the connection between 
galaxy interactions and star formation or AGN activity. Restricting 
the selection to pairs involving only early-types (i.e. precursors of 
dry mergers) results in a cleaner sample, minimising the contami- 
nation from gas in the interaction process. 

Within the standard framework of a ACDM cosmology, ellip- 
tical galaxies are formed through th e merger of many smaller sys- 
tems (see e.g. |Pe Lucia et al. 2005) and although th e build up of 
the red population occurs m ainly at high redshift (e.g. Bund v et al. I 
l2005l :l Ferreras et al. ]|200i), signs of previous merging events are 
found in early-type systems, such as kinemat ically decoupled cores 
jPavies et aUl200ll ; iMcDermid et al. Il2006h, distorted morpholo- 
gies, shells and other fine structure ( Ivan Dokkumll2005t) . 

The overall old stellar populations found in elliptical galax- 
ies can be used to constrain the time since the last major episode 



2 B. Rogers et al. 



10kpc | 10kpC 

• I 

■ b 



Figure 1. An example of close pairs of SDSS early-type galaxies usewd in this paper: Non-interacting (a; SDSS J 101434.27-00501 3.1; z=0.045) vs. interacting 
(b; SDSS J083645. 89+4722 10. 1 ; z=0.053). The vertical bar represents a projected distance of 10 kpc. 



of star formation. If the last merger involved late type galaxies, it 
must have taken place at early times, otherwise the resulting young 
populations would be seen in a photometric/spectroscopic analy- 
sis. However, if the last merger involves gas poor galaxies, one 
could still accommodate mergers at late times. It has been shown 
that models describing the formation of massive elliptical galax- 
ies require a nondissipative merger (i.e. limited amounts of cold 
gas) to recreate the dynamics of t he largest e llipticals, such as boxy 
isophotes and mi nimal rotation (Na ab et al. 1120061) . (although see 
iKang et al. l2007t) . Therefore it is reasonable to assume that the last 
merger for a significan t fraction of ellipticals inv olved early-type 
galaxies as progenitors (Khochfar & Burkert 2003). 

Given the large fraction of elliptical gal axies likely to con- 
tain sm all but significant amounts of cold gas ( Knap p"et al. II 19891 ; 
lYound l2005h . disruption of this gas via gravitational interac- 
tions/harassment or through 'dry' mergers may well result in small 
amou nts of star formation. It is also worth consid ering that previous 
work (Ferr eras et al. ^ 2006; Rog ers et al. 1120071) has indicated that 
recent star formation may be more significant in medium density 
environments, implying that galaxy harassment may be efficient at 
stimulating small episodes of star formation. 

In this paper we investigate these possibilities by looking at el- 
liptical only close pairs. We analyse the emission spectra to identify 
active galaxies (AGN or star forming) and explore whether the im- 
pending merger affects the nature of the activity. We also perform a 
comprehensive analysis of the stellar populations, using two inde- 
pendent methods: principal component analysis, and a maximum 
likelihood analysis involving a grid of models. The latter uses a 
number of age and metallicity sensitive spectral features measured 
via a new estimator of equivalent width that minimis es the contam- 
ination from neighbouring lines (Roge rs et all2008l) . 



2 THE CLOSE PAIRS SAMPLE 

Our classification starts with an initial sample of ~ 3000 galaxies 
in close pair s ystems from the Sloan Digital S ky Survey (SDSS) 
DR6 database ( Adelman-McCar thv et al. 120081) . A close pair is de- 
fined as a system of two or more galaxies within a projected phys- 
ical distance of 30 kpc and with relative velocity below 500 km/s 
(i.e. Az = 0.0017). The choice is motivated by previous work on 
the prese nce of interaction signs in close pairs of SDSS galaxies 
(see e.g. IPatton et all2000l) . 




0.05 0.1 1 10 



Redshift R p /R t 



Figure 2. The distribution of the sample is shown in terms of redshift (bot- 
tom left), central velocity dispersion, <tq (top left), apparent r-band magni- 
tude (top right) and the ratio between the projected separation of the pair 
and the half-light radius, Rp/R e (bottom right). The histograms show the 
distributions of the pairs depending on whether the system is visually clas- 
sified as interacting (black solid) or not (grey). The full sample corresponds 
to the black dashed lines. Overplotted are the cumulative distributions (with 
the same colour coding) which highlight differences between the two types 
of pairs. 



We make a visual classification of this sample using r-band 
FITS stamps retrieved from the SDSS pipeline. The morphology of 
each system of galaxies was estimated by eye by B. R., I. F., S. K. 
and A. P. In order to simplify the classification method, we bin the 
systems into three types: non-early type, lenticular or elliptical. We 
emphasize that our sample targets systems where both galaxies are 
early-types. Non-early type systems are those which show one or 
more galaxies with structure characteristic of a late-type: a domi- 
nant disk component or spiral arms. Lenticulars and ellipticals are 
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comparing the equivalent widths of absorption features with syn- 
thetic models in order to constrain the star formation history (SFH). 

Figure [2] shows the general properties of our early-type close 
pairs sample. The histograms are coded with respect to the interac- 
tion type: non-interacting systems are shown in grey, and interact- 
ing pairs are shown in black. The full sample is shown as a dashed 
histogram. We include in the figure the cumulative distributions to 
help assess the difference between samples. We draw attention to 
the central velocity dispersion (do; top-left), whose distribution is 
slightly biased towards higher values for interacting galaxies, as ex- 
pe cted {Park & Choil20 09). We note that a similar trend was found 
bv lAlonso et al. I ( 120070 ■ who showed that the estimated black hole 
mass, derived from the cto-BH mass relation, was systematically 
higher for interacting pairs. We find no difference with respect to 
interaction in the distribution of absolute magnitude (top-right) or 
redshift (bottom-left). While one could expect a bias in the classifi- 
cation, such that galaxies with brighter apparent magnitudes would 
have more visible debris, therefore increasing the fraction of "inter- 
acting" types, the distribution of apparent r-band magnitude shows 
no such trend (top-right). The distribution of interacting galaxies 
appears clearly biased towards smaller separations (defined as the 
ratio between pair separation and the half-light radius of the galaxy: 
Rp/R e ; bottom-right). 



Figure 3. The variation of emission line activity with respect to pair separa- 
tion is shown for our sample of close pairs (dashed lines, including Poisso- 
nian error bars) and for a control sample of SDSS early-type galaxies from 
ISchawinski et al. I i2007a. solid lines). For each panel, the left (right) side of 
the box corresponds to all galaxies (galaxies separated from the other mem- 
ber by less than 10 times the half-light radius). The filled (open) circles 
correspond to galaxies visually classified as interacting (non-interacting). 
This figure illustrates the fractions shown in tablefT] 



defined as bulge-dominated galaxies with or without a visible disk 
component, respectively. This classification resulted in a sample of 
346 early-type only pairs and one triplet system, comprising 695 
galaxies in total. 

These pairs were further divided into categories, based on 
the estimated level of interaction occuring within th e system. This 
schem e is similar to previous classifications (see lAlonso et al. I 
120071) . although we only consider two cases. Our systems are ei- 
ther "non interacting" or "interacting" depending on whether the 
close pair shows tidal tails or other distortions, indicating a signif- 
icant level of interaction. Figure [T] shows an example of each type. 
This criterion is also determined by a majority selection from the 
four classifiers. Systems for which there was an equal split between 
classifiers was considered as "non interacting". 

Optical spectra for the entire sample were retrieved from the 
SPSS DR6 and processed using the PPXF and G ANDALF codes 
( Cappel lari & Emsellemll2004l : Sarzi et al. 2006), adapted for fit- 
ting SDSS data. These codes provide stellar kinematics and emis- 
sion line lumino sities by simultaneou sly fitting both stellar pop- 
ulation models dTremonti et al. 1 120041) and gaussian emission line 
templates. The latter - when present - allow us to determine the na- 
ture of the ionising source. A comparison of the B aimer decrement 
(Hq/H /3) to standard case B recombination ( Ost erbrock & Ferlana 
2006) is used to correct for the effect of internal dust by using the 
extinction law of Calzetti et al. (2000). The emission spectra is re- 
moved from the observed data to determine the original SED of 
the underlying stellar populations. This is an important point when 



3 COMPARING THE EMISSION LINE SPECTRA 

Mergers and interactions disturb the gas in galaxies, and may lead 
to both AGN activity as the gas is driven towards the central super- 
massive black hole, or star formation if gas clouds can cool and col- 
lapse. Thus, identifying those galaxies which have undergone such 
activity in our sample will allow us to gauge the efficiency of this 
process in the specific case of gas-poor encounters between early- 
type galaxies. The identification of the ionising source in a galaxy 
is possible using emission line diagnostic diagrams (or BPT dia- 
grams, see Bald win et aL~lll98lh . In this paper we use the standard 
BPT diagram juxtaposing the [Nil] /Ha and [OIII]/H/3 line ratios 
dVeilleux & Oster brock 1987). The classification sc heme of galax- 
ies ac cording to this diagram is fully explained in iKewlev et al. I 
d2006h and is summarised below. 

Star-forming galaxies are found below the Kauffmann et al. 

d2003l) line: 

1Og([OIII]/H ^- lo g ( [ NII]/H Q )-0.05 +1 - 3 - (1) 

Galaxies above this line in the BPT diagram are likely to con- 
tain emission lines from an AGN. These galaxies are subdivided 
into two further sets: 'Composite' galaxies, possibly containing 
both an AGN and star forming regions, and galaxies whose ionisa- 
tion spectra is dominated by AGN activity. Composite objec ts are 
those galaxies which lie above the iKauffmann e t al. ( 2003]) star 
formation line but below the theoret ical upper limit of emission 
stimulated by star formation, namely (Kewl ev et al. Il200"lh : 

MPIII]/^) = log([NII] /Ha)-0.47 + ^ (2) 

Galaxies above this upper limit are said to have their emission 
dominated by an AGN. While it is possible to continue splitting the 
AGN population into Seyferts and LINERs, this adds little to our 
analysis, so it is omitted. Note that for a reliable classification, the 
four emission lines involved are required to have S/N 3. We iden- 
tify ~ 20% of our close pairs sample (~ 140 galaxies) as 'active' 
galaxies. 
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Table 1. Emission line classification for close pairs 



Sample 




Close Pairs 




Interacting 


Non interacting 


Comparison 1 




ALL 


R P /i? e < 10 


ALL 


Rp/fl e < 10 


ALL 


Rp/i?e < 10 




Number 


695 


140 


467 


41 


228 


99 


15729 


Quiescent 


83.4% 


79.3% 


82.0% 


76.8% 


84.1% 


85.4% 


81.5% 


Star Forming 


3.9% 


3.6% 


2.2% 


4.0% 


4.7% 


2.4% 


4.3% 


Composite 


3.2% 


1.4% 


3.1% 


2.0% 


3.2% 




6.9% 


AGN 


9.5% 


15.7% 


12.7% 


17.2% 


7.9% 


12.2% 


7.2% 



The comparison sample is taken from Schawinski et al. Id2007al) . 
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Figure 4. Top: The luminosity of the [OIII] (5007 A) line is shown with re- 
spect to the ratio between projected separation and half-light radius, Rp/R e 
(left), and stellar mass ratio (right) of the Seyfert, LINER and unclassified 
AGN galaxies. Galaxies in pairs classified as interacting (non-interacting) 
are shown as black solid (open grey) circles, respectively. The histograms 
are colour coded in the same way. Bottom: Specific star formation rates 
with respect to the ratio between projected separation and half-light radius 
(Rp/Re, left), and stellar mass ratio (right) for the star forming sub-sample. 
Galaxies in pairs classified as interacting (non-interacting) are shown as 
black solid (open grey) circles, respectively. The histogram is the distri- 
bution of spe cific star forma tion rates from the general sample of SDSS 
galaxies from lLi et al. 1 1200 8|) . 



3.1 Frequency of Active Galaxies 

Shown in table Q]are the fractions of the sample with signs of AGN 
activity as a function of the pair interaction type. For a control sam- 
ple of early-t ype galaxies we can not use our original set of early- 
type galaxies dRogers et al."l20 07~). That sample was extracted from 
iBernardi et al. I (2006), which incorporates a colour cut in the selec- 
tion, resulting in a bias on the number of active galaxies. Therefore 
we make use of a more general data s et comprising morphologi - 
cally classified early-type galaxies from lSchawinski et al. 1 12007a), 
taken from SDSS DR4 and spaninng a similar range in redshift. 
Figure [3] shows the dependence of emission line activity with re- 
spect to pair separation (Rp). The dashed lines illustrate the differ- 
ence between the full sample (left side of each panel) and a subsam- 



ple of close pairs, defined as those galaxies for which Rp < 10R e 
(right side of each panel). Short dashed lines show the Po issonian 
error bars. The control sample of lSchawinski et aT71 d2007al) is given 
by a horizontal solid line in each panel. We also illustrate the de- 
pendence on signs of interaction: filled (open) circles correspond to 
galaxies classified as interacting (non-interacting). One can see that 
the differences between early-type galaxies in close pairs and in the 
control sample is quite small, but it appears to be more significant 
in systems with visual signs of interaction. 

Previo us studies have found an increased level of both star 
formation dLambas et al. 1120031 ; lAlonso et al. 1120041 ; IWoods et al. 



| 2006t iKnapen & Jamesl 120091) and AGN activity dAlonso et al 
l2007UKeellll996l) within close pair samples not segregated with 
respect to morphology. However, table[T]shows that restricting the 
analysis to early-type galaxies results in little or no significant dif- 
ference - within Poissonian fluctuations - in activity between close 
pairs and the general population. This effect may be due to the lack 
of a substantial amount of gas in early-type galaxies to fuel either 
AGN activity or star formation. The small but non-negligible frac- 
tion of 'active' galaxies found both in a general and a close pairs 
sample could suggest that this activity is caused by a mechanism 
other than major mergers. A minor merger scenario dKavirai et al. I 
2009) could explain those fractions and the insensitivity to being in 
a close pair. 

On the other hand, while the emission line properties of 
the overall close pair sample is consistent with a general sam- 
ple, galaxies with visual signs of interaction show a h igher frac- 
tion of AGN compared to both lSchawinski et al. I d2007al) and our 
non-interacting subsample. This result is significant within Pois- 
son error b ars. We note t hat th e observed increase is in agreement 
with both Alons o et al. involving close pairs from SDSS, 

and iKavirail d2009). who targeted Luminous InfraRed Galaxies 
(LIRGs). The underlying fractions in both samples are clearly de- 
pendent on the type of galaxy targeted, but the increase is con- 
sistent. The other point to notice is that the close pair early-type 
galaxies with visual signs of interaction also feature a lack of star 
formation. This is in contrast to previous studies on close pairs of 
all morphological types, which find a strong link bet ween decreas- 
ing separation and increased star format ion (see e.g. lAlonso et al. I 
|2004| ; (Woods et al. ll2006l ; lLi et al. 1120081) . The interesting aspect is 
that the decrease in the fraction of star forming galaxies is coin- 
cident with an increase in the AGN fraction, a topic that will be 
discussed later on. Furthermore, the number of galaxies in the tran- 
sition region between AGN and star formation (bottom-right panel 
of figure [3} decreases at small pair separation (Rp < 10R e ). Al- 
though speculative, one could envision those trends as an increase 
in AGN activity with decreasing pair separation. 

Alternatively, it is possible that the signature of an increased 
black-hole accretion triggered by galactic interactions is damped by 
the presence of diffuse LINER-like emission that can masquerade 
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Figure 5. The luminosity of the [OIII] line as a function of velocity disper- 
sion. The sample is segregated into interacting (black) and noninteracting 
(grey). We note that the interacting sample has an increased average lumi- 
nosity but also that the correlation with mass is very weak for this sample. 



as AGN activity, as recently suggested bv lStasinska et al. 
This would be the case, in particular, if the power of such diffuse 
emission is limited by the amount of ionising photons produced 
by the old stellar populations that can reach the gas in early-type 
galaxies, rather than by the amount of such material. 

3.2 Star Formation vs. AGN activity 

To determine the possible effects of an increasing inflow of material 
into the centres of early-type galaxies in close pairs, we compare 
the star formation rate and the AGN activity from the luminos- 
ity of targeted emission lines. One would expect those to depend 
on the relative velocity, mass ratio or separation between galax- 
ies. The star formation rate (SFR) is estimated from the luminosity 
of the Ha emission line, using the standard correlation dKennicuttl 
1998). In order to factor out the stellar mass, one can also define the 
specific star formation rate (i.e. the SFR per unit st ellar mass) tak- 
ing th e stellar masses from the Garching catalogues jGallazzi et al. I 
120051) . The AGN activity can be traced using the luminosity of the 
[OIII] 5007A line which scales with the bolometric luminosity of 
the AGN an d thus the accretion ra te of the the central super-massive 
black hole (Heck man et al. |2004» . 

In figure[4](top) we show the luminosity of the [OIII] line as a 
function of the separation (Jeft) and the stellar mass ratio between 
the members of the pair (right). We only show those galaxies classi- 
fied as pure AGN, since massive stars will also affect the luminos- 
ity of the [OIII] line, contaminating the interpretation of L([OIII]). 
The accretion rate is found to correlate with separation. The corre- 
lation holds even if we consider only the non-interacting pairs (grey 
open circles), which are overall found at separations £ 10R e . This 
suggests that even in the absence of visual signs of interaction one 



Figure 6. Sample properties when segregated with respect to emission line 
activity. The solid black, dashed grey and dotted grey histograms corre- 
spond to quiescent, star-forming and AGN galaxies, respectively. 



could detect the effects of a close pair interaction from the activity 
of the central super- massive black hole. On the other hand, L[OIII] 
does not correlate with the stellar mass ratio (right). The vertical 
histogram (top-right) illustrates the trend between visual signs of 
interaction and the activity of the central nucleus. 

It is important to notice that this correlation does not reflect 
any bias regarding stellar mass. Shown in figure[5]is the luminosity 
of the [OIII] line as a function of velocity dispersion (i.e. a proxy 
for mass), we see that the increase of L[OIII] due to the mass of the 
galaxy is minimal compared to the correlation seen in figure|4]with 
respect to pair separation. 

As shown in table [TJ only 4% of the total sample have emis- 
sion lines consistent with star formation. Those galaxies are shown 
in the bottom panels of figure [4] where we present the specific star 
formation rate with respect to projected separation (Jeft) or mass 
ratio (right). On the far right panel, the distribution of the specific 
star formation r ate for the general sample of SDSS galaxies from 
lli et al. (2008) is given for comparison. As expected, our early- 
type close pairs sit at the low end of the distribution for a general 
sample. No significant trend is seen, although one could glimpse 
a correlation in the sample of non-interacting pairs (open grey cir- 
cles) such that the specific star formation decreases with decreas- 
ing separation. Although weak, this trend fits with the increasing 
L([OIII]) as the pairs get closer, suggesting a transition from on- 
going star formation towards AGN activity. 

Figure[6]shows the distribution of some properties of the sam- 
ple when segregating with respect to emission line activity. The 
solid black, dashed grey and dotted grey histograms correspond to 
quiescent, star forming and AGN galaxies, respectively. Given the 
large differences in the number of galaxies belonging to each class, 
we show there histograms normalized to the number of members in 
each class. It is worth noting the strong correlation with respect to 
galaxy mass (bottom-left), with the AGN and star forming galaxies 
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Figure 7. The distribution of the sample within r\ and f PCA space. The 
close pair sample (black histograms) is separated with respect to the visual 
classification: non-interacting (bottom) and interacting (top). For compari- 
son, SDSS early-type galaxies with a GALEX detection are shown as grey 
histograms, with the dashed (solid) lines corresponding to red (blue) galax- 
ies. The colour separation is based on the criterion for the presence of recent 
star formation, name ly NUV— r<4.9 for blue galaxies and NU V— r>5.9 for 
red galaxies (see e.g. lKavirai et a l. 2007; Rogers et al. l2007h . 



dominating the top and bottom ends in stellar mass, respectively. 
The characteristic mass of a quiescent galaxy sits in between these 
two. Regarding mass ratios (top-right), we find a significant en- 
hancement of star formation for more equal (but low-mass) ratios. 
This enhancement is also seen for lower relative velocities (defined 
as the difference between radial velocities; bottom-right). 



4 COMPARING ABSORPTION LINE SPECTRA 

In addition to AGN activity and ongoing star formation - which 
are processes triggered on shorter timescales - we also explore the 
effect of the interaction on the stellar populations as seen from the 
continuum and the absorption line spectra. In principle, one may 
not expect any significant changes in the properties of the bulk of 
the stellar populations during the first stages of an interaction (i.e. 
the phase we are only sensitive to in this sample). However, small 
and recent episodes of star formation can be detected in the opti- 
cal spectra as shown in otherwise red and dead early-type galaxies 
(Roge rTet al. l2007h . Galaxies with recent but no on-going star for- 
mation will be excluded from the analysis in the previous section, 
based on the emission line spectra. Hence, by extending our analy- 
sis to the continuum and the absorption lines, we increase our sen- 
sitivity to detecting the effect of a close encounter between early- 
type galaxies. We consider two independent methods to determine 
differences in the stellar populations. The first one involves Prin- 
cipal Component Analysis (PCA), a model-independent technique 
aimed at extracting "directions" in a vector space spanned by the 
spectral data, along which the variation is maximal. A more con- 



Figure 8. The parameters r] (bottom) and £ (top) derived from PCA are 
plotted as a function of the ratio between projected separation and half- 
light radius (Jeff); mass ratio (centre) and stellar mass (right). The black 
solid (grey open) circles are averages for the pairs with (without) visual 
signs of interaction. The vertical lines span the 25 th — 75 th percentiles in 
each bin. 



ventional second approach targets age-sensitive spectral features 
and compares them with a grid of models combined with popula- 
tion synfImsi£_STCCJra_to constrain the star formation histories (see 
e.g. lRogers et al. 1120081) . 



4.1 Principal Component Analysis (PCA) 

In order to maximally extract information from spectral data and 
to identify the smallest differences between them, we apply the 
method of Principal Component Analysis. It has proven to be a 
useful to ol in the analys i s of e xtre mely homogeneous sa mples, as 
shown in lFerreras et al. I d2006t) and lRogers et al. I d2007l) . where it 
was possible to identify small amounts of recent star formation at 
the level of a few percent in mass of ~1 Gyr old stars from the 
optical spectra. 

We take the principal components (i.e. the basis "spectral vec- 
tors" on which the SEDs of these galaxies are projected) from a 
general (i.e. non close-pairs) sample of ~7000 early-type galax- 
ies extracted from SDSS, previously defined and analysed in 
iRogers et al. I J2007h . We refer the interested reader to that paper 
for details of the method. This comparison set is a volume-limited 
sample ex tracted from the large r set of SDSS early-type galaxies 
defined in lBernardi et al. (2006), with the constraints: M r ^ —21 
and z^0. 1, with a further constraint on the signal-to-noise ratio of 
the spectra, S/N^ 15 per pixel. This control sample provides us 
with the pre-processed basis spectra for the analysis (i.e. the prin- 
cipal components). After de-redshifting and correcting for Galaxy 
dust absorption, the spectra from the early-type galaxies in close 
pairs are then projected onto these eigenvectors, giving the pro- 
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Figure 9. Best-fit stellar populations for a 2-burst model. The analysis involves a number of age- and metal-sensitive line strengths (see text for details). From 
top to bottom, the Mg-to-Fe ratio, metallicity, average age and the ratio between the young mass fraction and its age, are shown with respect to projected 
separation (left), mass ratio (centre) and central velocity dispersion (right). Solid black (open grey) circles correspond to galaxies with (without) visual signs 
of interaction. The vertical lines span the 25 th — 75 th percentiles of the distribution within each bin. 



jected components (i.e. PCI, PC2) that quantify the relative weight 
of each eigenvector in the construction of the spectra: 

N 

PCl i = $ i -ei=^$ i (Aj)ei(Aj) ) (3) 

3=1 

and similarly for PC2, PC3, etc... While it is possible to determine 
the principal components straight from the close pairs data set, the 
smaller sample size means the extracted eigenvectors will be less 
robust. PCA relies on the variance of the data set to define an eigen- 
vector and so large sets are preferred. 

4.1.1 rj and £ components 

The most successful mapping of PCA to extract information from 
the underlying stellar populations was found through th e first two 
principal components. As shown in lRogers et al. I d2007l) these two 
components, PCI and PC2, are consistent with an old and young 
stellar population, respectively. The eigenvector of PCI shows a 
pronounced 4000A break, significant metal absorption (e.g. H&K 
Ca II lines or a prominent G Band) and little or no Balmer absorp- 
tion. In contrast, the eigenvector of PC2 features a blue continuum 



with well-defined Balmer lin es. The positive cor relation between 
PCI and PC2 (see figure 3 of iRogers et al. 1120071) . is likely caused 
by the relative values required to reconstruct the spectra (e.g. the 
shape of the continuum), which depends on the age and metallicity 
of the stars. Relative to this relationship between PCI and PC2, a 
galaxy with a higher value of PC2 suggests the presence of a young 
sub-population. This idea was confirmed through a two-component 
stellar population model, in which the excess of the PC1-PC2 re- 
lationship was found to correlate with the mass fraction in young 
stars. 

We also found a consistent correlation between the projections 
of the principal components and NUV photometry from GALEX. 
The NUV spectral region (A ~ 2300A ) is very sensitive to the 
presence of small fractions of young stars. The (NUV— r) colour 
has been s hown to serve as an excellent indicator of r ecent star 
formation dSchawinski et al. l2007ahlKavirai et al. 1120071) . We used 
GALEX photometry to define two subsets of galaxies, the first one 
is NUV bright (NUV— r$;4.9) and represents galaxies that have 
undergone recent star formation (within ~lGyr). The second one 
is NUV fai nt (NUV— r^5.9) an d corresponds to an old, quiescent 
population. Rog ers et al. I d2007l) showed that the projections of the 
principal components - which only use the optical spectra - could 
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Figure 10. Colour-stellar mass relation of our sample. The colours are ex- 
tracted from the SDSS DR7 database, dereddened and measured over the 
3 arcsec aperture of the fibers used for the spectroscopy. The sample is 
divided with respect to visual signs of interaction (top panel) or emission 
line activity (bottom). In the top panel, solid black (grey open) circles cor- 
respond to galaxies with (without) visual signs of interaction. In the bot- 
tom panel, solid black (grey open) correspond to galaxies with AGN (Star- 
forming) activity. The grey crosses are quiescent galaxies. 



be used to discriminate between these two populations of galaxies. 
Thus from the linear fit of the PC 1 v. PC2 correlation, we define r\ 
as the distance along the relationship (sensitive to average age and 
metallicity) and £ as the residual from this fit (sensitive to recent 
star formation). 



4.1.2 Results 

Figure[7]shows the distribution of components r\ (left) and £ (right) 
for the close pairs (solid line histograms). The top (bottom) panels 
correspond to the interacting (non-interacting) visual type, respec- 
tively. In order to compare with a 'control' sample, we also show 
in each panel the histograms for the s ample of (non-close pair) el- 
liptical galaxies (Ro gers et al. 1120071) . segregated with respect to 
NUV— r colour. NUV bright (faint) galaxies are shown as grey solid 
(dashed) lines. Both "interacting" and "non-interacting" galaxies 
have the same distribution of the r\ and £ components. Thus, PCA 
indicates that there is little difference between the average popula- 
tions of the two classes, which suggests that the visual level of dis- 
ruption does not dictate the amount of recent star formation. Given 
these results it is not obvious that the interaction plays a signifi- 
cant role in shaping the stellar populations. This might be expected 
since our systems are the precursors of a dry merger. Furthermore, 
one could expect that the effect will not be apparent in the optical 
spectrum until later on in the merger. 

However, from the histograms on the right-hand side of fig- 
ure [7] one can see that the close pairs sample, both interacting and 
non-interacting, has a distribution of the £ component consistent 



with the general (i.e. non-close pairs) subsample of NUV-bright 
galaxies. This result indicates that elliptical galaxies in close pairs 
are more likely to contain small amounts of recent star forma- 
tion than a general sample of early-type galaxies. It is acknowl- 
e dged that this re s ult co uld be due to the selection effects present 
in Bernardi et al. 1 feOOfj) . However, the dearth of close pair galax- 
ies at the low end in the distribution of the £ component means 
that this bias is unlikely. The distribution in the r\ component is not 
so useful to discriminate between quiescent spectra and recent star 
formation, although the figure shows that the close pair distribu- 
tion sits 'roughly' in between the histograms for the NUV bright 
and faint galaxies. Given that 77 correlates with colour, one could 
conclude that studies purely based on broadband photometry can- 
not measure these differences. We note that these distributions are 
in agreement with previous work on early-type galaxies in Hickson 
Compact Groups, where a larger scatter towa rds high values of PC2 
(roughly Q were found JFerreras et al. 2006). 

In figure [8] we show the PCA components 77 and £ as a func- 
tion of pair separation (Jeft); stellar mass ratio (centre) and stellar 
mass (right). The black solid (grey open) circles correspond to the 
interacting (non-interacting) pairs, and the vertical lines span the 
25 th — 75 th percentile of the distribution in each bin. Out of the 
three observables, stellar mass is the only one that correlates sig- 
nificantly with the components, appearing redder (i.e. more neg- 
ative 77) towards higher masses. The correlation is much weaker 
with respect to f , which is the proxy for the presence of recent star 
formation. The mass ratio does not present any significant trend, 
and only the non-interacting pairs give a slight trend towards an 
increased recent star formation at small separations, although the 
scatter of the sample is rather large. 



4.2 Modelling the Star Formation History 

Since the close pairs sample has a higher mean value of the £ 
component, we expect a large fraction of the sample to have un- 
dergone recent star formation. In order to quantify the effect on 
the underlying stellar populations, we explore a two-component 
star formation history. The synthetic spectra are generated from 
the 2007 stellar population models of Ch ariot & Bru zual (see 
e.g.lBruzual & Charlotll2003l ; |Bruzualll2007h . assuming a lChabrieJ 
(2003) initial mass function. The composite model superposes two 
simple stellar populations: an old component with age to, al- 
lowed to vary between 2 and 14 Gyr, with the metallicity between 
log(Z/Z@) = —1.5 and +0.4. A younger component of the same 
metallicity is added, with age tY, contributing a mass fraction, fy, 
which ranges from to 0.5. The age of the young component is 
taken between 100 Myr and 2 Gyr. Note that the model grid con- 
tains a subset of models which are equivalent to standard simple 
stellar populations (i.e. fy = 0), such that a composite model will be 
chosen only if it improves the fit. Also note that the oldest ages con- 
sidered (14 Gyr) are motivated by our choice of a standard ACDM 
cosmology (Sl m =0.3, Ho = 70km/s/Mpc). The final grid of model 

consists of 65,536 star formation histories. 

Our analysis follows the approach of iRogers et al. I d2008l) . 
which involves multiple age-sensitive spectral features comprising 
three Balmer lines (H/3, H7 and H8) and the 4000A break strength 
(D4000), along with a metal-sensitive index, [MgFe] (as defined in 
iGonzalezll 19931) . Furthermore, a new definition of equivalent width 
(EW) is used, that significantly reduces the age-meta llicity degen- 
eracy over the traditional side-band method (see e.g. iTrager et al. I 
2000). Our EWs are based on a new definition of the pseudo- 
continuum, determined from a (boosted) median of the surrounding 
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spectra and a 20A spectral window centered on the line of inter- 
est. This definition has been shown to reduce the contamination of 
the pseudo-continuum from neighbouring lines, resulting in a less 
metal dependent H7 and H<5 and a less age dependent [MgFe]. This 
method also provide s smaller uncertainties in the EW at low S/N 
dRogers et all2008l) . 

While the BC03 models are not c alibrated to ac c omm odate 
non-solar abundance ratios, following Yamada et al. ] J2007h wc 
use as a proxy the comparison between the best fit values of the 
metallicity when replacing in the analysis the metal-sensitive index 
[MgFe] by either (Fe) - defined as 0.5(Fe5270+Fe5335) - or by the 
Mg6 index. The difference in the metallicity derived from these two 
fits - labelled as Ziig and Zf c - is given as a crude estimation of the 
abundance ratio: [ZMg/ZFc] = log(ZMg/ZF C ). The comparison of 
observed data and models is done via a standard maximum like- 
lihood method. The errors for the EWs are estimated from Monte 
Carlo realizations of gaussian noise applied to the spectra. These 
errors are adde d in quadrature to the es timated systematic errors 
for the models dBruzual & Charlotll2003l) . 

4.2.1 Results 

The results of the line strength modelling are shown in figure [9] as 
a function of projected separation (Jeft); stellar mass ratio (centre) 
and central velocity dispersion (right). Black solid (grey open) cir- 
cles correspond to close pair galaxies classified as interacting (non- 
interacting). The error bars give the 25 th — 75 th percentile range 
of the distribution within each bin. From top to bottom, we show 
the proxy for abundance ratio ([ZMg/Zp a ]); metallicity; (mass- 
weighted) average age, and a quantity that gives the strength of the 
recent episode of star formation. The mass and age of the young 
population is degenerate, such that a small mass in young stars can 
be replicated by a larger mass fraction in older stars. Hence, we 
parameterise the effect on the spectra in terms of the ratio between 
the mass fraction in the young component to its age (Ty/Iy). 

Consistently to the PCA studies described above, there is no 
significant difference between the stellar populations with respect 
to the visual presence of interactions. Furthermore, the popula- 
tions are only sensitive to the velocity dispersion, a result already 
prese nt in general samples of galaxies (see e.g. iBernardi et al. I 
There is no significant trend with respect to separation or 
mass ratio. However, the modelling of the line strengths do reveal in 
an independent and consistent way to PCA that early-type galaxies 
in close pairs are more likely to have undergone recent star forma- 
tion. We find 378 out of 695 galaxies (i.e. 54%) have a significant 
young population, log(fY/tY) ^ —2, which is equivalent to 1% 
mass fraction of a 1 Gyr population. 

The correlation between colour and stellar mass can also be 
used to understand the connection between pair morphology, emis- 
sion line activity and the underlying stellar populations. Figure [Tol 
shows the (u — r) colour-stellar mass diagram, where the colours 
are directly obtained from the SDSS DR7 database (de-reddened 
and measured within the 3" aperture of the spectrograph fibers). In 
the top panel, solid black (grey open) circles correspond to galax- 
ies with (without) visual signs of interaction. In the bottom panel, 
solid black (grey open) circles correspond to galaxies with AGN 
(star forming) activity. The grey crosses are quiescent galaxies. As 
expected, colour is strongly correlated with mass, which confirms 
the trend seen in the PCA components in figure [8] However, the 
departure from the red sequence - which is an indicator of younger 
stellar populations - does not depend on signs of visual interaction, 
agreeing with the PCA result of figure [7] However, with respect to 
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Figure 11. The distribution of the sample is shown in terms of estimated 
mass of the dark matter halo of the occupied galaxy group (Mh). The his- 
tograms show the distributions of the pairs depending on whether the sys- 
tem is visually classified as interacting (black solid) or not (grey). The full 
sample corresponds to the black dashed lines. 

emission line activity (bottom panel of figure [Tot, we find a signif- 
icant trend such that quiescent galaxies and AGN systems (crosses 
and solid dots, respectively) populate the red sequence, whereas 
galaxies with on-going star formation are identified as the mem- 
bers of the 'blue cloud'. The evolution from a weak star forming 
cloud towards a weak AGN and subs equently a quiescent gal axy 
has already been proposed elsewhere ( Scha winski et al. l2007ah . In 
the conclusion we apply this evolutionary path to the interpretation 
of the early-type close pairs. 



5 ENVIRONMENT 

We can also explore the properties of the close pairs with respect 
to environment. In order to assess the environment of a galaxy, we 
use the estimated mass of the host dark matte r halo. This is deter - 
mined through the group finding algorithm o f lYang et al. I d2005l) . 
which identifies galaxy groups starting with a friends-of-friends al- 
gorithm. The membership to these groups follows an iterative pro- 
cess controlled by the properties of the group and its halo. The 
application of this algorithm to a larg e sample of SDSS galaxies, 
forms the galaxy groups catalogue of lYang et al. I d2007l) . The in- 
terested reader is directed to those references for an in depth de- 
scription of the algorithm and catalogue. A cross correlation of this 
catalogue to our close pairs sample reveals a match to 664 of the 
695 galaxies. As a measure of environment, we use the halo mass 
estimated based on the observed stellar mass of the group. We note 
that a more detailed investigation into the effects of environment as 
measured by the mass of the galaxy group halo is also underway on 
a much larger sample of ellipticals (Rogers et al. in preparation). 

FigureQj]shows the distribution of our sample in terms of the 
host halo mass. The galaxies from our close pairs sample are mainly 
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located in intermediate mass halos of M_s ~ IO^Mq/i -1 . This is 
in agreement with the halo mass found by IPasquali et al. I d2009h 
at which galaxies change from being preferentially star-forming to 
showing optical-AGN activity. 

Our sample is clearly too small to map global properties of 
active early type galaxies, however it is interesting to see how 
the activity of t he clo se pairs are affected. In a manner similar to 
IPasquali et all {2009), we show in figure [T2] the conditional frac- 
tions of galaxies classified as having undergone significant recent 
star formation (top) or split with respect to their star forming / AGN 
activity (bottom). The conditional fractions are given by the num- 
ber of active galaxies within a bin of host halo mass divided by the 
total number of galaxies in that bin. The dashed lines towards the 
left of the figure represent the halos for which a considerable bias is 
expected towards low mass galaxies as imposed by the mass of the 
host halo. In addition to the limitation of galaxy mass by the group 
mass itself, it should be reme mbered that low m ass halos preferen- 
tially host low mass galaxies dYang et al. 120081) . In a sample of this 
size it is difficult to overcome such a bias and this caveat should be 
considered alongside the conclusions of this section. 

Figure [T2l (bottom panel) shows that the fraction of galaxies 
classified as AGN is fairly constanQup toM ff ~ 4xlO 13 M /i _1 , 
above which there appears to be a rapid decline. The dearth of AGN 
galaxies at high h alo m as ses is consisten t with p revious results 
iKauffmann et al. I fe004l) ; iGilmour et ail d2007l) : IPasquali et ail 
J2009h . The exact reason for the cut off is not obvious but may 
be related to the reduction of gas available due to increased tidal 
stripping. 

The fraction of star forming galaxies is higher in lower mass 
halos, although low mass g alaxies - well known to have higher 
fractions of star-formation ( Kauffma nn et al. I l2003h - dominate 
within these bins as mentioned above. The rest of the sample show 
a consistent fraction of star forming galaxies across all environment 
types. The number of galaxies is relatively small so a robust con- 
clusion cannot be inferred, but it is interesting to note that no sig- 
nificant drop in star forming galaxies is seen with respect to halo 
mass. 

Also shown in figure[T2](iop panel), is the effect on the stellar 
populations where we consider how the amount of recent star for- 
mation is affected by the mass of the host halo. We define a "signif- 
icant" amount of recent star formation as Sy = log(f Y /t Y ) ^ —2, 
or £ ^ 10 -3 . The fraction of galaxies with RSF drops by ~10 - 
20% as we move from low/intermediate to high halo masses. 
We note that this drop i s qualitatively similar, alt hough relatively 
higher, than that found in lSchawinski et al. I ( I2007bl) . which may be 
due to the high sensitivity of NUV or a function of the enhanced 
RSF already seen in the sample. 



6 CONCLUSIONS 

We have selected a sample of ~350 close pairs involving only 
early-type galaxies to measure environment effects. Through the 
use of emission line diagnostics, we classify ~ 20% of the sam- 
ple as containing either an AGN or currently undergoing (weak) 
star formation. The fraction of these 'active galaxies' is consis- 
tent with an independent sample of early-type systems not in close 



1 The small or zero fraction of AGN at low halo masses is most likely due 
to the lack of intermediate and high mass galaxies in these halos. However 
the increase in SF in these bins may hint that dust from such activity may 
obscure some of the AGN. 
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Figure 12. (Top Panel) The change in the fraction of galaxies with signifi- 
cant recent star formation as a function of the host halo mass (Mg). Shown 
for both PCA and SFH modelling results, where the significant level of RSF 
is defined by C-3 s C/10' 3 ^ 1 and (Sy = ) log(f Y /t Y ) ^ -2. (Bot- 
tom Panel) The fraction of galaxies containing AGN or star formation as a 
function of environment. The fraction are conditional, f(AGN, SF\Mn), 
refering only to the galaxies contained within the specific halo mass bin. 
The solid grey and black histograms correspond to star-forming and AGN 
galaxies, respectively. The dashed line signify bins containing a bias to- 
wards low mass galaxies. It is interesting to note the drop in AGN fraction 
after log M H ~ 13.6 (M /i -1 ). 



pairs, taken from ISchawinski et al. I d2007al) . However, we found 
an excess of AGN and a lack of star forming galaxies in close pairs 
with visible signs of interaction. This result suggests that during the 
encounter, galaxies evolve from a (weak) star forming phase to an 
AGN phase. This idea is supported by the increase both in the AGN 
fraction and in the luminosity of the [OIII] (5007A) line towards 
decreasing pair separation. Additionally, the specific star formation 
rate shows a hint of a decrease with decreasing separation, although 
the number of star forming systems is very low (contributing only 
4% to the total sample). 

A significant increase is found in the number of early-type 
galaxies in close pairs that have undergone a recent star forma- 
tion episode, with respect to a control sample. This is shown both 
through PCA, where the close pair sample features high values of 
the £ component (sensitive to young stellar populations), as well 
as by the line strength analysis using a grid of 2-component mod- 
els, where a large proportion of the sample requires a significant 
amount of young stars. These two apparently contradicting scenar- 
ios - namely that the encounter appears to reduce ongoing star for- 
mation and the observation that the sample has increased levels of 
recent star formation - can be reconciled. 

It is speculated here that the observed recent star formation is 
triggered during the first phases of the encounter by the interactions 
of the outer parts of the galaxies, including the dark matter halos. 
Since a large fraction of ellipticals have been found to contain HI 
and molecular gas, not only in their interstellar medium, but also in 
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the form of satellite gas clouds jKnapp etal. II 1989c Young 2005; 



iMorganti et al. ll200dTcombes et al. ll2007l ; lDonovan et al. II2007I) . 
the increasing gravitational perturbations induced by an oncom- 
ing neighb our will destabilise these clouds, driving them to wards 
the galax y dSofue & Wakama7svJll993l ; ldi Matteo et al. 1120071) . For 
instance, iLi et al. I d2008h found an enhanced star formation rate 
up to separations of 100 kpc on a large sample of SDSS galaxy 
pairs. Indeed the accretion of gas has been s hown to instigate star 
formation even in early - type galaxies (e.g. ISofue & WakamatsvJ 
ll993l;|Pipino et al. pOOa I Khalatvan et al. 1120081) . The simulations 



of Khalatv an et al. I ( 120080 suggest that the accretion of gas onto 



an early-type galaxy will be short lived due to the feedback from 
AGN. 

As the pair comes closer - and within our selection criterion of 
Rp < 30 kpc - the increased gravitational interaction will enable 
the rem oval of angular momen tum, driving gas to wards the centre 
(see e .g. iD'Ercole et al. ||2000|) . Simulations from di Matte o~et al. I 
1 20071) indicate that the greatest inflow of material to the centre 
occurs at separations ~10 kpc. Hence, at lower separations, most of 
the available gas is driven towards the centre, triggering the AGN 
activity and possibly q uenching star formation . This is consistent 
with the simulations of Ijohansson et al. I J2008h . in which mergers 
between elliptical galaxies showed decreasing star formation rates 
with progession of the merger and rapid termination at the later 
stages coincident with increased black hole accretion. The feedback 
from the AGN in such cases should drive out the majority of the 
gas within the galaxy into the intergalactic medium. Given the high 
metallicity of early type stellar populations we might expect this to 
contribute to the abundance of the IGM. 

Although speculative, this scenario explains the transition 
seen in the emission and absorption features of our sample, and 
is consistent with the observations of lSchawinski et all j2007ah on 
a general sample of early-type galax ies. Our results agre e with the 
general study of SDSS close pairs of Elliso n et al. I d2008l) who find 
that star formation also precedes AGN activity in close pairs of 
late-type galaxies. 

The environment in which the interaction takes place also 
seems to affect the exact nature of the encounter. In low mass ha- 
los the initiated star formation appear to continue for longer into 
the interaction, i.e. to appear within our selection window of 30 
kpc separation. The amount of recent star formation also appears 
to be higher in these low mass halos (possibly as a consequence) 
as well as in intermediate mass halos, which harbour most of the 
AGN activity. In contrast, halo masses greater than Mjj ^ 4 x 10 13 
Mq/i -1 appear to have less AGN activity (almost none) and re- 
duced amounts of RSF. 
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